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ABSTRACT 
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Context. BL Lac objects show core-jet structures with features moving outwards along the jet. We present a kinematic analysis of jet 
component motion in the pc-scale jet of the BL Lac object S5 1803-1-784, which does not reveal long-term outward motion for most 
w- \ . of the components. 

S»i/ ■ Aims. S5 1803-F784 shows complex kinematic phenomena; understanding these provides new insights into the emission processes in 

BL Lac objects and possibly into the differences between quasars and BL Lac objects. 

Methods. The blazar S5 1803-1-784 has been studied with VLBI at v=1.6, 2.3, 5, 8.4, and 15 GHz between 1993.88 and 2005.68 in 26 
observing runs. We (re)analyzed the data and present Gaussian model-fits. We collected the already published kinematic information 
Q ■ for this source from the literature and re-identified the components according to the new scenario presented in this paper Altogether, 

^-j ' 94 epochs of observations have been investigated. 

[/3 I Results. A careful study of the long-term kinematics reveals a new picture for jet component motion in S5 1803-1-784. In contrast to 

, previously discussed motion scenarios, we find that the jet structure within 12 mas of the core can most easily be described by the 

coexistence of several bright jet features that remain on the long-term at roughly constant core separations (in addition to the already 
known "stationary" jet component ~ 1.4 mas) and one faint component moving with an apparent superluminal speed (~19c, based 
on 3 epochs). While most of the components maintain long-term roughly constant distances from the core, we observe significant, 
smooth changes in their position angles. We report on an evolution of the whole jet ridge line with time over the almost 12 years 
of observations. The width of the jet changes periodically with a period of ~8-9 years. We find a correlation between changes in 
the position angle and maxima in the total flux-density light-curves. We present evidence for a geometric origin of the observed 
phenomena and discuss possible models. 

Conclusions. We find evidence for a significantly different scenario of jet component motion in S5 1803-F784 compared to the 
generally accepted one of outwardly moving jet components, and conclude that the observed phenomena (evolution of the jet ridge 
' line, roughly constant component core separations but with significant position angle changes) can most easily be explained within a 

' geometric model. 

Key words. Techniques: interferometric - BL Lacertae objects: individual: S5 1803-1-784 - Radio continuum: galaxies 

. 1 . Introduction However, as more detailed VLBI-monitoring reveals, this model 

d ' can not easily explain all the observed motion phenomena in a 

■ ■ ■ . Jet component motion m most Active Galactic Nuclei (AGN) ^^^^^^^ ^^^^^^ ^qN. In many sources stationary features 

has been successfully described by relativistic shock-m-jet ^^^^ ^^^^^^ ^^^^^^ ^^e core have been observed (e.g., 

motion directed away from the presumed central engine (e.g., ^^^^^^^ 200O; Kellermann et al. 2004; Britzen et al. 2007). 
Marscher & Gear 1985). Bulk relativistic motion explains many 
observed jet-related phenomena and is also important with 

regai-d to the observed short-term radio variability (IDV, Wagner Taking the core as the easternmost component, the brightest 

& Witzel 1995 and references therein) providing the necessai-y jet feature (in 8.4 GHz observations) at r = 1.4 mas used 

Doppler factors to reduce the observed brightness temperatures to be one of the most prominent candidates for a stationary 

to the Compton Hmit (Kellermann & PauUny-Toth 1969). component (e.g., Eckart et al. 1986; Schahnski et al. 1988). 

However, as we showed in an earlier work, this so-called 



Send ojfprint requests to: Silke Britzen Stationary component has a varying core-separation (Britzen 

* Member of the International Max Planck Research School (IMPRS) et al. 2005a). Based on a more detailed analysis of new data, 

for Radio and Infrared Astronomy at the Universities of Bonn and we show in this paper that most of the jet components in the 

Cologne pc-scale jet of S5 1803H-784 (out to r ~ 14 mas) seem to remain 
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at roughly constant core separations on the long term, but their 
position angles undergo significant change. We present the 
results of a detailed kinematic analysis of the pc-scale jet in this 
source and discuss possible explanations of these phenomena. 



1.1. TheblazarSS 1803+784 

As a member of the complete S5 sample of 13 flat-spectrum 
radio sources at high declinations (Witzel 1987), S5 1803H-784 
has been observed repeatedly at many frequencies and angular 
resolutions since the late 1970s (e.g., Eckart et al. 1986, 1987; 
Schalinski et al. 1988; Witzel et al. 1988; Krichbaum et al. 
1993; Britzen & Krichbaum 1995). These observations revealed 
the complex morphology of the milli-arcsecond-jet of this 
intra-day variable source (Wagner & Witzel 1995 and references 
therein). The source shows a pronounced jet with prominent 
jet components located at core separations of r = 1.4, 5, and 
12 mas (e.g., Eckart et al. 1986). Geodetic and astronomical 
VLBI data gathered at 8.4 and 5 GHz between 1979-1987 
showed that the component located at r - 1.4 mas appears 
stationary (e.g., Schalinski et al. 1988, Witzel et al. 1988). 
Several authors confirm this constant core separation, leading 
to a subluminal velocity of ySapp < 0.74/;"' (e-g-, Cawthorne et 
al. 1993). Britzen et al. (2005a) presented an overview over the 
pc-scale structure of S5 1803H-784, including discussion of the 
kinematics of this stationary r = 1.4 mas component. Based on 
a jet-component identification procedure in which the brightest 
jet component in every epoch is associated with this one, this 
r = 1.4 mas component approaches the core whenever a new jet 
component was or was about to be ejected from the core. Thus 
this stationary component was found to have non-constant core 
separation ("oscillatory" type behaviour). Britzen et al. (2005a) 
explain this motion scenario within a model of reconfinement 
shocks. In addition, three jet components were reported to move 
with apparent superluminal speeds of 7- 10c. 

S5 1803H-784 shows curved jet structure on all scales 
probed with VLBI. The curvature found in the cm-regime is 
described in Britzen et al. (2005b). In the mm-regime, the inner 
jet components move superluminally with expansion rates of 
0.14 + 0.04 mas yr"' and 0.07 ± 0.05 mas yr"', with possibly 
variable proper motions (0.02-0.28 mas yr"^) along a curved 
path, suggesting helical motion (e.g., Krichbaum 1990). This 
source also belongs to the group of misaligned AGN (Antonucci 
et al. 1986; Appl et al. 1996), with the pc- and kpc-scale jet 
aligned almost perpendicularly. The probable transition region 
between the pc- and kpc-scale jet has been investigated in world 
array observations by Britzen et al. (2005b). A summary of the 
observations of the kpc-scale structure can be found therein. 

Velocities computed in the previously cited papers have 
been determined with differing sets of cosmological parameters. 
Here, we use the following: Ho = 71 km s^'Mpc"', Q.m= 0.27, 
Qa= 0.73, and take z = 0.68 (Stickel et al. 1993; Lawrence et al. 
1987) 

1.2. New scenario for jet motion in S5 1803+784 

This is the third paper in a series of papers on the radio structure 
of S5 1803H-784 and presents an alternative explanation for the 
observed kinematics in the pc-scale jet compared to the scenario 
presented in Britzen et al. (2005a). Preliminary results have 



already been published in Kudryavtseva et al. (2006). 

An obvious difference between the motion scenario pre- 
sented in this paper compared to that in Britzen et al. (2005a) 
is the significantly larger time span covered by the observations 
upon which the present one is based. By combining observa- 
tions covering more than 20 years (in the case of the 1.6 GHz 
observations), it is possible to investigate the long-term motion. 
This long-term behaviour shows no significant outward motion 
of most of the components in this source. 

The most prominent jet component in the pc-scale jet of S5 
1803H-784 at frequencies between 5 and 22 GHz is the 1.4mas 
component to the west of the easternmost, brightest component. 
This r = 1 .4 mas component has repeatedly been characterized 
as stationary. As densely sampled geodetic observations have 
shown, this stationary component approaches the core from 
time to time, showing a sort of oscillatory behavior (Britzen 
et al. 2005a). In this paper we take note of the fact that since 
Eckart et al. (1986) other jet components have repeatedly been 
observed at similar core separations of e.g., 6 and 12 mas 
(Britzen et al. 2005b and references therein). 

Despite the fact that some jet components have been 
observed at roughly constant core separations over long periods, 
they have been identified in the literature with components 
showing long-term outward motion. We here propose a new 
component identification assuming that the components ap- 
pearing repeatedly at the same distance from the core indeed 
represent the same component. We compare our results with 
results published in the literature and apply our identification 
scheme to these data and find that the literature data are fully 
consistent with our new identification scenario. 

Several stationary components have been reported for a 
number of other AGN in the literature. These components re- 
main at a similar position in core separation and position angle. 
However, here we encounter jet components that remain at 
roughly constant core separations but show significant changes 
in their position angle. 

We use multi-frequency data to study the kinematics of the 
components, present a correlation analysis among a number of 
properties including the flux-density evolution, propose an al- 
ternative identification scenario, and conclude with the result- 
ing implications. The structure of this paper is organized as fol- 
lows: section 12] overviews the data sets that have been observed 
and analyzed or re-analyzed. Section[3]presents the results of a 
detailed model-fitting analysis of the 26 epochs of observation 
and shows evidence for an alternative identification scenario and 
unusual jet kinematics in S5 1803H-784 based on 94 epochs al- 
together (data re-analyzed plus literature data). In addition, we 
compare our results with kinematic information derived from the 
literature. Section|4]summarizes the kinematic results. Section 5 
presents correlations and periodicities from our analysis. Finally, 
section|6]briefly discusses the implications for our understanding 
of kinematics and emission processes in AGN. 

2. Observations and data reduction 

We have (re-)analyzed 26 multi-frequency (v=1.6, 2.3, 5, 8.4, 
and 15 GHz) VLBI observations of S5 1803H-784, performed 
between 1993.88 and 2001.87 with heterogeneous arrays: VLBI 
observations at 1.6 and 2.3 GHz conducted by Marcaide et al. 
(1995a, 1995b, 1997); at 5 GHz by Marcaide et al. (1995a, 
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1995b, 1997), Guirado et al. (2001), and Gurvits et al. (priv. 
comm.); at 8.4 GHz by Perez-Torres et al. (2000), Ros et al. 
(2000, 2001, in prep.); and at 15 GHz by Perez-Torres et al. 
(2000, in prep), Kellermann et al. (1998), Zensus et al. (2002), 
Lister & Homan 2005. The data have been fringe-fitted and 
calibrated before by the individual observers. For calibration 
details we refer the reader to the original publications. 

In addition, we performed space- VLBI observations at 1.6 
GHz in 2001.87. These observations were scheduled for the 
VLBA plus HALCA. Unfortunately we did not obtain any 
data from ground-space baselines and use this observation 
as a purely ground-based experiment. The data have been 
fringe-fitted within AlPS. 

We fit circular Gaussian components to each of the data sets 
at each frequency using the difinap package (v.2.4b. Shepherd 
1997). In order to find the optimum set of components and pa- 
rameters, we fit every data set starting from a point-like model. 
We used circular Gaussian components in order to avoid the ef- 
fects of correlations among weakly-constrained axial-ratio es- 
timates and those of other parameters of the Gaussian compo- 
nents. 

2.1. Determination oftfie uncertainties 

The determination of model-fit parameter uncertainties is a com- 
plex topic and no satisfying solution to this problem has been 
presented so far. A detailed discussion of model-fitting uncer- 
tainties and their influence on the parameters is presented in 
Britzen et al. (2007, 2008). The uncertainties of each data set 
in general are influenced by different u-v coverages, calibra- 
tions, sensitivities, etc. The relative positional errors are different 
at different core separations and the brightness of the individ- 
ual component affects the fitted position as well. Uncertainties 
should reflect all these different error sources. In order to ob- 
tain reliable error estimates we used two different algorithms to 
determine the model-fitting uncertainties. In the first, the un- 
certainties for the core separation and the position angle are 
given by the formulas Ar = (c/cr„„ , -y/1 + S peaki crrms)l'2.S peak and 
A0 = arctan(Ar/r), where cr„„5 is the residual noise of the map 
after the subtraction of the model, d is the full width at half max- 
imum (FWHM) of the component and S peak is flie peak flux den- 
sity (Fomalont 1999). However, this formula tends to underesti- 
mate the uncertainty if the peak flux-density is very high or the 
width of the component is small. In the second, we calculated all 
uncertainties by comparing different model fits (+1 component) 
obtained for the same set of data. These uncertainties reflect the 
possible parameter ranges for the individual components within 
model-fitting. To be conservative, we selected the maximum of 
the uncertainties from the two methods as the value of the uncer- 
tainty to use in subsequent analysis. 

2.2. Data table 

In Table [T] we list the parameters of the best fits to the data and 
their uncertainties. We also give the component identification 
and the reference to the original publication of the data sets. In 
Figs. [T5l420l we show the hybrid images with model-fit compo- 
nents superimposed. For the images from observations between 
February 2001 and September 2005, the beam sizes have been 
adjusted for better comparison. The original beams were: epoch 
2002 (15 GHz) 0.4207 x 0.6068 mas, epoch 2003 (15 GHz) 



0.4764 X 0.4929 mas, epoch 2005/2 (15 GHz) 0.4768 x 0.4962 
mas, epoch 2001.09 (8 GHz) 0.697 x 1.036 mas, epoch 2001.29 
(5 GHz) 1.137 x 1.651 mas. 

2.3. Component identification 

Our kinematic analysis relies critically on the component iden- 
tification. Our new scenario is based on the flux-density, core 
separation, and position angle of the individual components. All 
components have been identified in a way such that all of these 
parameters vary as little as possible between adjacent epochs. 
For those components that have been identified at several fre- 
quencies, we did not find a significant frequency-dependent off- 
set. The results are presented and discussed in detail in the fol- 
lowing sections. 



3. Results 

3. 1. Lacl< of long-term r(t) cfiange in most of ttie model jet 
components 

In addition to the formerly so-called stationary component (Ca) 
at ~ 1.4 mas, we find components at core separations of ~0.3 mas 
(CO), ~ 0.8 mas (CI), ~2 mas (C2), ~3^ mas (C4), -6-8 mas 
(C8), and ~10-12 mas (C12). Most of these jet components — 
based on the new data presented here — remain at similar core 
separations between 1993.88 and 2005.68. The core separations 
of C12 and C8 at 15 GHz are less well defined, reflected by their 
larger uncertainties. Since the identification of C12 suffers from 
a reduced number of data points, we do not include this com- 
ponent in our further analysis. Components C30 and C40 can 
only be traced at 1.6 GHz, and also will not be considered in the 
detailed analysis. Fig.[T]panels (a-d) present an overview of the 
r(f) behavior for all jet components at 15, 8.4, 5, and 1.6 GHz, 
respectively, out to a core separation of ~50 mas and over a time 
span of up to twenty years. Panels (b) and (c) show data col- 
lected from the literature, but with component re-identifications 
following the scenario we present in this paper; these previously 
published data support this new component identification sce- 
nario. 

3.1 .1 . Comparison with data from the literature 

In order to check whether the components we find showed 
a similar or different kinematic behavior in earlier times, we 
coflected afl available VLBI information for S5 1803H-784 from 
the literature. We list the model-fit parameters derived for the 
total intensity VLBI observations in Table |2] We re-identified 
all the components from these published model fits according 
to the scenario proposed in this paper and plot the results in 
Fig. [T] panels (b) and (c). At both frequencies (5, 8.4 GHz), 
the components have remained at similar core separations 
for almost 20 years; we find no long-term outward, apparent 
superluminal motion based on the aggregate data. Although 
different kinematic scenarios have been proposed by different 
authors, the data are fully consistent with a scenario in which 
the components tend to remain at long-term roughly constant 
core separations. 

Panel (d) of Fig. [T] shows r(t) from the 1.6 GHz data. 
We again find that the inner components (up to 12 mas) 
maintain roughly constant core separations, although the outer 
components (C30 and C40) do show some evidence for outward 



4 



S. Britzen et al.: The kinematics in the pc-scale jets of AGN 




Fig. 1. Panel (a) shows the core separation as function of time, r(f), for those jet components detected at 15 GHz. Individual 
components are denoted by different symbols and lines (usage common across all panels). Note the straight line for component Ca 
and the fast moving component B3. Panels (b) and (c) show r{t) for 8.4 GHz and 5 GHz observations, respectively, based on data 
collected from the literature, with component re-identification following the new scenario discussed in this paper Panel (d) shows 
r(f) for components farther out along the jet based on 1.6 GHz observations (new + from the literature). Panel (e) shows the motion 
of component Ca in rectangular coordinates at 8.4 and 15 GHz, in which a displacement towards larger x values is apparent between 
the fits from the earlier (1986-1994) and newer (1995-2001) 8.4 GHz data with the 15 GHz fits falling in between. Panel (f) shows 
the flux-density along the jet based on all new epochs of the 15 GHz data. The x and y axes scales differ across panels (a)-(d). 



motion. 

In Table |3] we give the parameters of the linear regres- 
sions performed for the jet components' core separation as a 
function of time, r^ean and fir and calculated apparent speeds. 



ySapp. Based on the values for the apparent speeds calculated 
for all frequencies, we find no superluminal motion for the 
components within ~4 mas of the core. CI and C2 show no 
significant motion within the errors. CO and Ca show some 
evidence for subluminal outward motion. Although C4 shows 
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apparent superluminal motion at 15GHz, combining data at 
all frequencies gives an apparent speed consistent with zero. 
We find significant apparent superluminal outward motion of 
the faint component B3. C8 and C12 show some evidence 
of apparent inward motion. For these two components fewer 
epochs having well-constrained positions are available. In the 
following we concentrate on the inner components. 



3.1 .2. Comparison with geodetic VLBI data 

Britzen et al. (2005a) presented an analysis of the kinematics in 
the pc-scale jet of S5 1803H-784 based on geodetic VLBI data. 
Such geodetic observations are performed in general more fre- 
quently than astronomical observations (please find details in 
Britzen et al. 2005). In order to have comparable numbers of 
observations per time from astronomical and geodetic observa- 
tions, we smoothed the geodetic model-fit results at 8.4 GHz 
in the period 1986.21-1993.95 with the following sliding mean 
algorithm. The new data points were produced, replacing the 
set of circular Gaussian parameters m,- for each time f, with the 
weighted mean: 

1 * 

— Ypjmj, (1) 

where k is the number of data points within interval - A, -i- A] , 
and the weight of jth point is determined as 

Pj^exp[-i6tj/Af], (2) 

where Stj is the time span from the jth point to the center of 
the window. For our analysis we selected a window value of 
A - 0.3 years, which is consistent with the temporal resolution 
of data points obtained at other frequencies. 



3.1 .3. Flux-density evolution along the jet 

Fig. [T] panel (f) shows the 15 -GHz flux-density as function of 
distance along the jet out to r = 2 mas, with the various epochs 
of observations overplotted. We find maxima in the observed dis- 
tribution at ~0.9 mas and 1.4 mas. For the component at 1.4 mas 
we find a trend towards a sharp peak in flux-density. 

3.1 .4. Displacement of Ca 

Fig. [Upanel (e) shows the position of component Ca in rectan- 
gular coordinates. The different symbols denote frequency (8.4 
and 15 GHz), with the 8.4 GHz results further divided into pre- 
and post-1995.0 sets. The latter 8.4 GHz data are obtained from 
the literature. These two 8.4 GHz data-sets seem to have a shift 
in right ascension of ~0.3 mas. We are comparing here geode- 
tic data with astronomical data. In the case of the geodetic data 
(Britzen et al. 2005a), we assumed that the brightest component 
is identical with component Ca. We investigate this effect in the 
scenario of a Binary Black Hole in more detail (Roland et al. 
2008). 

3.2. The "fast" component m 

In addition to the stationary components, one fast moving 
component (B3) has been observed in the time between 2002.1 




-? -d -e -8 



Fig. 2. Model fits superimposed on the hybrid maps obtained 
from Cleaning and self-calibration for the most recently ob- 
served five epochs. The component positions are connected by 
dotted lines. The positions of the fast component B3 are con- 
nected by a solid line. Please note that the elapsed time between 
adjacent epochs differs. 

and 2005.7. As shown in Fig. [T] panel (a) and Fig. |2] this 
component shows outward motion with an apparent velocity 
of ~19c based on 3 epochs of observations. B3 is much fainter 
(between 0.01 and 0.04 Jy) and shows a kinematical behavior 
that differs significantly from the other components, which 
don't show long-term secular outward motion. The ejection 
time of this component computed by extrapolating it's proper 
motion back to r = is 1999.8 ± 1.1. Unfortunately, this epoch 
is in the middle of the two-year gap in the observations. 



3.3. Further evidence for the coexistence of one "fast" 
component among "stationary" components 

In order to check our model of the S5 1803H-784 jet, in which 
only one component moves with apparent superluminal speed, 
we investigated in detail the opacity and spectral evolution of 
the total flux-density light-curves. The emergence of a new jet 
component, according to the shock-in-jet model (e.g., Marscher 
1996; Gomez et al. 1997), is caused by a shock induced by a 
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Table 3. Linear regression fits to the motion of listed components, r^nemi is the mean core separation over the time-range of the 
regression, /u is the proper motion of the component, and /Japp is the apparent speed, fi and /3app are computed twice: once using only 
15 GHz results, and once with the combined results from all frequencies, to is the extrapolated time of component ejection from the 
core in the case of component B3. 



Comp. ID 


'■mean [mas] 


fil^ [mas yr '] 


Alpp [f] 


fif [mas yr ' ] 




to [yr] 


CO 


0.30 ± 0.02 


0.017 ±0.001 


0.39 ± 0.02 


0.006 ± 0.001 


0.14 ±0.02 


/ 


CI 


0.79 ± 0.03 


-0.037 ± 0.003 


-0.86 ± 0.07 


0.001 ±0.001 


0.02 ± 0.02 


/ 


Ca 


1.27 ±0.03 


0.007 ± 0.004 


0.16 ±0.09 


0.022 ± 0.001 


0.51 ±0.02 


/ 


C2 


1.93 ±0.03 


0.017 ±0.007 


0.39 ±0.16 


-0.002 ± 0.001 


-0.05 ± 0.02 


/ 


C4 


3.70 ± 0.08 


0.145 ±0.017 


3.37 ± 0.39 


0.001 ± 0.001 


0.02 ± 0.02 


/ 


C8 


6.78 ±0.18 


-0.182 ±0.005 


-4.23 ±0.12 


-0.182 ±0.005 


-4.23 ±0.12 


/ 


C12 


10.02 ±0.38 


-0.715 ±0.113 


-16.62 ±2.63 


-0.197 ±0.069 


-4.58 ± 1.60 


/ 


B3 


/ 


0.807 ±0.151 


18.76 ±3.51 


/ 


/ 


1999.8 ± 1.1 



primary excitation at the base of the jet. This manifests itself in 
radio light-curves as an outburst that is delayed at lower frequen- 
cies due to the combined effects of the frequency stratification of 
the emitting electrons, non-zero opacity, and light-travel delays. 
Such time-delayed outbursts are associated observationally on 
mas-scales with the brightening of the VLBI core and are ac- 
companied by the ejection of a new optically thin jet component. 



We decomposed the total flux-density light-curves at 
4.8 GHz, 8 GHz, and 14.5 GHz from the University of Michigan 
Radio Observatory (UMRAO) monitoring campaign (Aller et 
al. 1999, 2003) into Gaussian components, as was described in 
Pyatunina et al. (2006, 2007). Fig. |3] shows the light-curve of 
S5 1803H-784 at 14.5 GHz and the resulting Gaussian decompo- 
sition. Three prominent outbursts are visible: (A) in 1985; (B) 
in 1988, and (C), a very prolonged outburst starts in ~1992 and 
ends in ~2005 with a peak in 1997. The Gaussian parameters, fit- 
ted to the light-curve are listed in Table |4] Frequency-dependent 
time delays for each flare are the time difference between the 
Gaussian peak at each frequency with respect to the position of 
the peak at 14.5 GHz. These are plotted in Fig.|4] Flares A and B 
show moderate time delays of 0.3 years, whereas the prolonged 
flare C shows an enormous time delay of 3.2 years between 
4.8 GHz and 14.5 GHz, which is an indication of high opac- 
ity in the source during the C outburst and a flat spectrum that 
becomes gradually steep. Outburst C has a very broad profile 
at all frequencies (0 ~ 8 yr). Similarly broad outbursts were ob- 
served for the source 1308H-326 in Pyatunina et al. (2007), where 
it was mentioned that such a broad outburst can be associated 
with the development of a dense region, some kind of "cocoon" 
that confines the innermost part of the jet that is excited by the 
primary perturbation. We calculated quasi-simultaneous spectra 
using the measurements at all three frequencies for which time 
separation is less than two weeks; Fig.|5]shows the spectral evo- 
lution fit as a power law to each spectrum. The spectral behav- 
ior changes significantly after ~1996, where the spectra change 
from being flat, and gradually become steeper. The changes of 
the spectral evolution occur in the beginning of a bright, pro- 
longed flare C and can be associated with it. We suggest that this 
can be explained by a jet component ejection. Small frequency- 
dependent time delays and steep spectra of other observed flares 
of the source suggest that there probably were no jet compo- 
nent ejections before 1995 which is in good agreement with our 
model for the S5 1803H-784 jet. 
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Fig. 3. Gaussian decomposition of the UMRAO 14.5 GHz fight- 
curve. The solid line shows the sum of the Gaussian functions 
and dotted lines are individual Gaussian functions. Letters indi- 
cate the names of the flares. 

3.4. Significant position angle clianges 

As shown before, most of the components remain at similar core 
separations for the time covered by the observations. Although it 
is obvious from the panels in Fig.[T](a) and (b) that most jet com- 
ponents show an additional motion with regard to the core sepa- 
ration superimposed on their rather stationary behavior, we con- 
centrate in this paper on the most important aspects only. A de- 
tailed analysis of the second order motion is complicated by the 
uneven sampling. This effect of "oscillating" motion has been 
analyzed in detail for Ca in 18034784 by Britzen et al. (2005a). 
The phenomenon of "oscillations" is clearly seen in other com- 
ponents of 18034784 and in other BL Lac objects as well and 
will be analyzed in detail in Glueck et al. (in prep). 
The components show significant position angle changes. Fig.|6] 
panels (a)-(d) show the evolution of the position angle for four 
jet components (CO, CI, C2, and Ca, respectively) at different 
frequencies. The position angles vary over 35° in the case of CO 
and 15° in the case of C2. We find very similar relations at 5 
and 8 GHz for CO, at 5, 8, and 15 GHz for CI, at 8 and 15 GHz 
for C2, and at 5, 8, and 15 GHz for Ca. One outlier for CO in 
Fig.|6]panel (a) in 1994 can probably be explained by a blending 
effect with the core or a new component due to the closeness of 
the core. In the case of the core separation for CO and C2, the 
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Comp. 


Freq. 


Amplitude 


Tmax 









Time delay 




GHz 


Jy 


yr 




yr 




yr 


A 


8.0 


1.94 ±0.02 


1985.19 


±0.02 


1.41 


±0.02 


0.00 ± 0.02 


A 


4.8 


0.77 ± 0.04 


1985.39 


±0.04 


1.41 


±0.04 


0.20 ± 0.06 


B 


14.5 


2.01 ± 0.01 


1988.93 


±0.01 


2.63 


±0.01 


0.00 ±0.01 


B 


8.0 


2.15 ±0.02 


1988.86 


±0.01 


4.07 


±0.02 


-0.07 ± 0.02 


B 


4.8 


1.25 ±0.01 


1989.19 


±0.01 


2.54 


±0.01 


0.26 ± 0.02 


C 


14.5 


1.46 ± 0.01 


1996.94 


±0.02 


9.83 


±0.03 


0.00 ± 0.02 


C 


8.0 


1.57 ± 0.02 


1998.34 


±0.04 


7.89 


±0.04 


1.40 ±0.06 


C 


4.8 


0.93 ± 0.01 


2000.13 


±0.04 


6.28 


±0.05 


3.20 ± 0.06 



3.5 




-0.5 ' ' ' ' ' ' ' 

4 6 8 10 12 14 16 

Frequency [GHz] 

Fig. 4. Time lags of individual outbursts as functions of fre- 
quency. 
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Fig. 5. Time evolution of the quasi-simultaneous spectral index, 
fit to the 4.8 GHz, 8 GHz, and 14.5 GHz flux densities. 



evolution is similar at the two frequencies. For C2 we find some 
indication that the lower frequency-data can be found at larger 
core separations. 



3.5. Evolution of the mean jet ridge line 

In Fig. [8] we show the jet ridge line for all epochs obtained at 
15 GHz in Cartesian coordinates. The jet ridge line is defined by 
a line connecting all component positions at a given epoch and is 
independent of any component identification. We find that the jet 
ridge line changes its shape between epochs in a periodic way. 
Fig.[8]panel (a) shows the temporal evolution of jet component 
positions at 15 GHz for the time period 1994-2005. It overplots 
the jet ridge lines at each epoch, with each dot representing the 
position of one jet component and each line connecting the jet 
components for one particular epoch. Panels (b) and (c) show 
this evolution with time in more detail; the left-hand columns 
showing the trace of the ridge line at each epoch, and the right- 
hand columns the flux density of the components along the jet 
as a function of x. From an almost straight line in 1994.67, the 
shape of the jet evolves into a sinusoid with the northernmost 
value at X ~ 1 mas and the southernmost value atx~2 mas. The 
amplitude of the sinusoid reaches its maximal values at epoch 
1998.84 and decreases again, forming an almost straight line 
again in 2003.10. One period of this excursion is completed af- 
ter ~ 8.5 years and the jet shape starts to evolve into a sinusoid 
again. However, the position of the straight lines in 1994.67 and 
2003.10 are different, with an offset in j of 0.1 mas. We thus 
see an evolution of the jet ridge line with time as defined by the 
positions of the individual components within the jet. 

3.6. Evolution of the jet width 

In this section we investigate the width of the jet as defined by the 
position angle distribution that the jets components span. Fig.|9] 
shows this distribution as function of time, including all jet com- 
ponents seen at all observed frequencies. This plot shows that 
the jet changes its width with time: around 1995 the jet width 
appears to be quite small (~5°), around 1999.5 the jet width 
appears to be much broader (~ 40°) and is narrower again in 
2004. To investigate the longterm evolution of this effect, we in- 
clude data taken from the literature in Fig. |9] panel (b) and find 
a similar evolution of the width. The complete dataset (our data 
+ literature data) seems to cover two cycles of jet-width broad- 
ening with a period of about 8-9 years. To investigate whether 
this broadening is mainly a frequency dependent effect, we show 
the distribution at 15 GHz in panel (c) and at 8.4 GHz in panel 
(d); both panels confirming the trend of changing jet width. To 
investigate whether less well-constrained model components at 
large core separations with large position angles produce this ef- 
fect, we checked whether the shape of the plot changes when we 
include only jet components within a certain range of core sep- 
arations. We checked the data for the jet components within the 
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inner 1 mas, 2 mas, etc. up to 7 mas. However, the shape of the 
position angle distribution is stable and does not change signifi- 
cantly. Panels (c) and (d) discussed above show the distribution 
for components within 5 mas of the core. 

To check whether periods of more frequent observations pro- 
duce more outliers in position angle and thus affect the position 
angle distribution, we binned the data in 0.3 year bins and cal- 
culated the standard deviation in position angle within each bin, 
both for 5 GHz and 8.4 GHz. This procedure used the compo- 
nents within 5 mas of the core. Fig. [T] shows the evolution of 
the standard deviation per bin (the outliers were removed from 
the plot), using the center of each bin as its epoch for plotting. 
The position angle spread with time shows a smooth increase 
starting in 1985, and reaches its maximum in 1988-1992. It then 
decreases again and has the second maximum in about 1998. 
It is clearly seen that the spread of the position angle evolu- 
tion has two cycles with a characteristic timescale of about 8- 
9 years, similar to the period suggested by the evolution of the 
jet ridge line. The autocorrelation function gives the value of the 
timescale: 9.9+0.2 with a correlation coefficient of 0.79. Binning 
the data with different bin length, such as 0.5 or 1.0 years gives 
the same results. In order to check whether the different num- 
ber of jet components for different epochs of observations can 
affect these results, we selected only data points within 0.5 mas 
core separation. We only took the first three jet components for 
each particular epoch into account and binned the position an- 
gles into 0.3 year bins. The evolution of the standard deviation 
of the position angles in these bins looks similar. We thus con- 
clude that the shape of the plot and characteristic time scale of 
the position angle changes do not depend on the number of the 
jet components for a particular epoch of observation. 
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Fig. 6. Position angle as function of time at different frequencies 
for four individual components. 
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Fig. 7. Time evolution of standard deviation of position angles, 
calculated for components within 5 mas of the core in 0.3 year 
bins. Dotted circles show 8.4 GHz data; circles 15 GHz data. 
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Fig. 8. The jet ridge line for all components detected in 15 GHz observations within 6 mas of the core. Panel (a) overplots all epochs, 
as denoted by the different lines/points. Panels (b) and (c) show the trace of the jet in Cartesian coordinates (left-hand columns) and 
the flux density along the jet as a function of x (right-hand columns). 
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4. Summary: kinematic resuits 

Before we describe the periodicities and correlations found 
in S5 1803+784, we give a brief summary of the results of 
the kinematic analysis. In Britzen et al. (2005a) we presented 
a component identification scenario for the pc-scale jet in 
S5 1803+784 based on the assumption that the brightest 
component seen in each observing epoch can be identified with 
the same component (Ca). The most plausible identification 
scenario resulting from this assumption was the one presented 
in that paper: the coexistence of several superluminally moving 
and one oscillating component. 

The data presented in this paper cover a broader range in 
frequency (1.6-15 GHz and time (1993.88-2005.68). Thus, 
a more detailed investigation of the jet component motion is 
possible. In addition, we collected kinematic information from 
the literature and tested our hypothesis against these data. In 
total, 94 data sets have been investigated. Conclusions derived 
from the present analysis include: 

- most jet components within the inner ~ 10 mas of the core 
remain at roughly constant core separation over the long term 
(see Fig.[T]i. 

- the position angles of these components change significantly 
over time (see Fig.|6]). 

- component (B3) moves with a superluminal apparent speed 
of ~ 19c based on 3 epochs (see Figs.[T]panel (a) and|2]i. 

- the jet ridge line evolves with time, with a likely period of 
~8.5 years (see Fig.|8]l, independent of any component iden- 
tification. 

- the width of the jet changes periodically with a similar 
period of ~8-9 years (see Fig.|9|l- 

In the following section we show that the epochs of the max- 
imum widths of the position-angle distribution correlate well 
with maxima in the total flux density observed within the 
UMRAO monitoring programme (see Figs. [TOl and fTTI). We find 
convincing correlations between parameters such as the flux- 
density, core separation, position angle, and describe them be- 
low. 

5. Correlations & Periodicities 

5. 1. Correlations & Anti correlations 

As shown in previous sections, the jet width changes with time 
and the jet ridge line itself shows a quasi-periodic variation 
with a characteristic time scale of ~ 8.5 yrs. Both facts could 
result from a geometric origin. This would also produce cor- 
relation among the core separation, the position angle, and the 
flux density changes for each jet component and across the 
different components. To investigate these correlations, we first 
performed a visual analysis of the time series of the position 
angle, core separation, and the flux of each component at 8.4 
and 15 GHz. Maxima and minima at 8.4 GHz are listed in Table 
|5] Major changes in the variations of the jet parameters at 8 GHz 
occurred in 1987, 1988.5, 1990 and 1991 for CO, CI, Ca and the 
core. The core flux reached a maximum in 1987, accompanied 
by minima in flux, position angle, and core separation of CI, 
maxima of the core separation of CO and maxima for flux and 
position angle of Ca. In 1988.5, the values that were maximal 
became minimal and vice versa, and switched again in 1990 



and 1991. This indicates that changes in flux, position angle 
and core separation of the components and the core-flux are 
correlated and change with characteristic time scales of a few 
years. Similar behavior was found as well at 15 GHz at two 
other epochs (Table |6l). The major changes in core separation, 
position angle, and flux variations of the inner components 
occun-ed in 1998 and 2000. 

To investigate the correlations among these circular 
Gaussian parameters of a component and across multiple 
components more carefully, we calculated the discrete cross- 
correlation functions (Edelson & Krolik 1988) for the pairs of 
parameters (core separation, flux density, position angle) for 
each of the components (CO, CI, Ca, C2, C4) at 8.4 GHz. 
Component C2 at 15 GHz possibly comprises 3 components at 
8.4 GHz. We cafl these three C2 components at 8.4 GHz C2i, 
C22, and C23. We find the following results: 



- For the brightest jet components Ca and CI, all the param- 
eters are correlated. As an example. Fig. [12] shows the time 
series of the core separation, the position angle, and the flux 
density of the component CI from 1984 to 1996. 

- For the jet components C2i, C22 and C4, we find a strong 
correlation or anti-correlation among the variations in flux 
density, position angle, and core separation within each com- 
ponent. Table|7]shows the discrete correlation function peaks 
for these parameters in jet components. 

- The core flux is correlated with the core separation, the flux 
density, and the position angle changes of all inner jet com- 
ponents CO, CI, Ca, C2i, C22 and C23 (see Table^i. 

- Variations of individual circular Gaussian parameters also 
show correlations across different jet components. Table |8] 
also shows the results of the discrete correlation function cal- 
culation. 



When cross-correlating core separation changes of Ca and 
those of other components, significant peaks in the discrete 
correlation function (e.g., CO, CI, C2i and C23) are associated 
with time lags that change gradually with distance along the 
jet. The time lags appear to increase roughly linearly with core 
separation; Fig. [14] shows the time lags for these component 
pairs as a function of the components' mean core separa- 
tion (time lags for components interior to Ca have had their 
sign inverted from Table |8]l, and fit with a linear regression. 
Similar behavior is found for the cross-correlation between 
the core separation changes of CI and Ca, C2i, C22, and 
C23. This can be explained if we assume that the motion of 
the components is mostly due to a geometrical origin. The 
components follow almost the same trajectories, but with a 
certain time delay. The components are connected with each 
other and the information is spreading with a proper motion 
of about 2.3 mas/yr if the reference point is component Ca 
and 0.79 mas/yr if the reference point is component CI, corre- 
sponding to extremely fast speeds of 53c and 18.2c, respectively. 

We investigated possible correlations among the jet- 
component parameters and the total flux-density variability. For 
this purpose we used the UMRAO data of multi-frequency AGN 
monitoring (Aller et al. 1999, 2003), spanning almost 30 years 
(see Fig. ^l^. Using the 8 GHz UMRAO light-curve over the 
time-range spanned by the VLBI data, we find: 
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Fig. 9. Position angle distribution of individual jet components as a function of time. Panels (a) and (b) show the new observations 
from this paper and those combined with the observations from the literature, respectively, for all frequencies. Panels (c) and (d) 
show components within 5 mas of the core 15 GHz and 8.4 GHz, respectively. 

Table 5. Maxima and minima in flux density, position angle and core separation at epochs 1987, 1988.5, 1990, and 1991 at 8.4 GHz. 



Time 


1987 


1988.5 


1990 


1991 


Comp. ID 


Flux 


P.A. 


Core Sep. 


Flux 


P.A. 


Core Sep. 


Flux 


P.A. 


Core Sep. 


Flux 


PA. 


Core Sep. 


Core 


max 






min 






max 






min 






CO 






max 






min 






max 


min 




min 


CI 


min 


min 


min 


max 


max 


max 


min 


min 


min 


max 


max 


max 


Ca 


max 


max 


min 


min 


min 


max 


max 


max 


min 









- variations in the core separation and position angle of Ca and 
C2i are correlated with the total flux-density variability (see 
TableEl). 

- variation in the core separation and the position angle 
changes of the jet component CI are correlated with the total 
flux density at 14.5 GHz (see Fig.[T2li. but do not show any 
correlation at 8 GHz. 

In order to investigate the same correlations at 15 GHz 
where the number of data points was not sufficient for the 
discrete correlation function analysis, we calculated the 
Pearson's correlation coefficient and inspected them visually. 
Tables Q and [8] show the Pearson coiTelation coefficients and 
the probabilities for finding such a value by chance for various 
pairings of different parameters for a single component and 
the same parameter across two components, respectively. The 
flux-density variability is anti-correlated for the pairs of Cl/Ca, 



Table 6. Maxima and minima in flux density, position angle and 
core separation at epochs 1998 and 2000 at 15 GHz. 



Time 


1998 


2000 


Comp. ID 


Flux 


PA. 


Core Sep. 


Flux 


PA. 


Core Sep. 


Core 


min 






max 






CO 


max 


max 


min 








CI 








max 


max 


max 


Ca 


max 


min 


max 


min 


max 


min 



and Ca/C2. We only find correlated core separation changes for 
the pair Ca/C2. Fig.[T3]shows notable examples for a coiTelation 
and anti-correlation: the position angle vs. core separation in 
component Ca (panel a) and the flux densities in components 
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Ca vs. CI (panel b), respectively. 
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Fig. 10. The spread in the position angle for four different fre- 
quencies as a function of time, based on all the data (this paper 
+ literature). We clearly find several maxima that are correlated 
with the maxima found in the total flux-density data. 



1803+784 
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8 GHz+ 1.5 Jy 
14.5GHz + 4Jy 




Time [yr] 

Fig. 11. Total single-dish flux-density light-curve for S5 
1803H-784 at 4.8, 8 and 14.5 GHz. The 8 and 14.5 GHz light- 
curves are shifted by 1 .5 and 4 Jy, respectively. All the data are 
from the UMRAO monitoring campaign. 



5.2. Periodicities 

To test for possible periodic changes in the component param- 
eter variations, we applied a discrete autocorrelation function 
(Edelson & Krolik 1988), a data-compensated discrete Fourier 
transform (Ferraz-Mello 1981), and the Jurkevich method 
(Jurkevich 1971). All three methods were applied over the time- 
range 1981-2005 for all components out to C12. We find a 
quasi-periodicity of about 2 and 4 years in the variation of the 
core separation, the position angle, and the flux with time for 




1984 1986 1 



1990 1992 1994 1996 
Time [yr] 



Fig. 12. Core separation, flux density, position angle of compo- 
nent CI, and total flux-density light-curve at 8 GHz in the period 
1984-1996. The two dashed lines indicate the nearly simultane- 
ous peaks in all four plots. 

the inner jet components CO, CI, Ca, and C2 at 8 and 15 GHz. 
This agrees well with the periods of 2 and 3.9 years found by 
Kelly et al. (2003) in the total flux density variability. A more 
detailed description of this analysis will be presented by N.A. 
Kudryavtseva et al. (in prep.). 

6. Discussion 

In this paper we have presented multi-frequency kinematic 
analysis for S5 1803H-784 based on observations covering 
roughly twenty years in time. The most important result is that 
most of the features in the pc-scale jet up to a core distance of 
~12 mas remain at roughly constant core separations over long 
time scales. This is in contrast to most previously presented 
kinematic scenarios for this source or flat-spectrum AGN 
in general. While many authors find that component Ca is 
stationary, most authors claim outward motion for the remaining 
components. In contrast, the 2cm/]V[OJAVE webpage and Lister 
et al. (2009) show a core separation/time plot in which several 
components within ~ 8 mas of the core also maintain roughly 
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Fig. 14. Time lag of the discrete cross-correlation function of the 
core separation of component Ca vs. those of components CO, 
CI, C2i, and C23 as a function of the mean value of distance 
from the core for these components. The Unear regression fit to 
these time lags is overplotted. 



separation, and the position angle are obviously related. Many 
other pairs of correlations or anti-correlations are presented in 
this paper. The set of all these correlations and periodicities to- 
gether argues in favor of a geometric origin. In the following sub- 
sections, we discuss possible scenarios to explain the observed 
evolution in the pc-scale jet of S5 1803H-784 and the kinematic 
behavior of its components. 



Fig. 13. Panel (a) shows core separation as a function of position 
angle for the brightest component Ca. Panel (b) shows the flux 
density of component CI against that of component Ca. Both 
panels plot the parameters estimated from 15 GHz observations 
from all the different epochs. 



constant core separations. This is consistent with our results 
presented in this paper. The 2cm data from the 2cm/M0JAVE 
survey form part of the hterature data sets included in the 
present analysis. The fast-moving component B3 is not seen by 
the 2cm/MOJAVE. 

The combination of new and re-analyzed data spanning 
12-20 years shows that all bright components of the inner jet 
remain at roughly constant core separations, with no significant 
outward superluminal motion. Applying our component identi- 
fication procedure to available data sets from the literature leads 
to a consistent kinematic picture. Although these components 
remain at roughly constant core separations on the long term, 
they do exhibit significant changes in their position angles. This 
results in an evolution of the jet ridge line from a straight line to 
a sinusoidal form with a characteristic time-scale of 8-9 years. 
A similar evolution of the jet ridge line has so far only been 
observed in one other BL Lac object: 0716H-714 (Britzen et 
al. in prep.). Following such jet ridge-line evolution proceeds 
independently of any component-identification-scenario but 
requires a large number of epochs to be discovered and traced 
over several periods of evolution. 

We find a number of correlations between the parameters of 
the core and those of individual components. The best case is 
seen for CI where the total flux, the component flux, the core 



6. 1 . Curved jet structures / Precession ? 

Curved jet structures have been observed in many different 
sources, both extragalactic (e.g., 1928H-738, OJ 287, 3C 273, 3C 
449) and Galactic (e.g. SS433, HHlll). There are prominent 
differences between stellar and AGN jets, including the role 
of magnetic fields in the collimation and dynamical evolution 
of their structure and the jet-to-environment density ratio that 
is significantly larger for the stellar jets. Based on the later 
characteristic, the dynamical evolution of a stellar jet can more 
easily be treated ballistically, whereas such an assumption 
would not necessarily apply to the case of AGN jets (Hardee et 
al. 1994). 

Curving and wiggling of jets is usually attributed to preces- 
sion of the base of the jet. Two different scenarios concerning 
the central engine that feeds the jets can explain the observed 
curvature. 



6.1 .1 . Models relying on a single central Black Hole 

In the first scenario, the core that ejects the jet is considered to 
be a single object (usually a supermassive black hole, SMBH). 
Precession can then occur due to the Bardeen-Petterson ef- 
fect through Lense-Thirring precession of the accretion disk 
(Caproni et al. 2006). However, this model assumes a rigid 
body configuration between the accretion disk and the jet, an 
assumption that we are not able to test or verify. A magnetic 
torque mechanism has also been proposed (Lai 2003), driven by 
the large-scale magnetic field threading through the accretion 
disk, leading to precession of the disk and consequently of the 
jet. In the same frame, a group of models attributes the curving 
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and wiggling of the jet to magnetohydrodynamic (MHD) insta- 
bilities, most prominently Kelvin-Helmholtz (K-H) instabilities 
(Camenzind & Krockenberger 1992; Hardee & Norman 1988; 
Birkinshaw 1991; Zhao et al. 1992; Hardee et al. 1994; Hardee 
et al. 1997; Meier & Nakamura 2006; Perucho et al. 2006). 
Meier & Nakamura (2006) consider the instability to be driven 
by current, related to the strong-field pinch, rather than by K-H 
instabilities. All models mentioned above strongly depend on 
the jet-to-environment density ratio. However, instabilities have 
to be excited in order to affect the structure of the jet. Thus, 
although the MHD instabilities (and mainly the helical K-H 
kink instabihty) can possibly explain the structure of pc and 
kpc scale jets, the problem of how these modes are excited stiU 
remains. 



6.1 .2. Models relying on a Binary Black Hole 

In the second scenario, models have been built under the 
assumption that a binary system can be found in the centers 
of AGN. The most popular idea concerns binary black hole 
systems (BBH) as a product of galaxy merging. The problem 
is complex and each of the models makes a number of as- 
sumptions to explain the observed structures. In such systems, 
the primary black hole is considered to have an accretion disk 
and jets, whereas the secondary black hole is in orbit around 
the primary. Gravitational perturbations associated with this 
orbiting secondary black hole can disrupt the accretion disk 
of the primary black hole and lead to precession. Lobanov & 
Roland (2005) argue that in light of the "two-fluid" model, the 
orbital motion and the disk precession can lead to curving and 
wiggling of the jet. Furthermore, they propose that for strong 
enough magnetic fields, the K-H instabilities can be considered 
negligible. Roos et al. (1993) argue that the orbital motion 
of the binary is responsible for the change of the jet-ejection 
direction and thus leads to a wiggUng structure (the jet is treated 
balUstically). Kaastra & Roos (1991) argue that the precession 
angle is determined by the angle between the primary black hole 
spin axis and the angular momentum of the binary system. As 
before, the jet is treated ballistically. Similar models have been 
proposed by Katz (1997) (a gravitational torque induced by the 
secondary black hole onto the accretion disk) and Romero et al. 
(2000) (near-rigid body configuration). 

Roland et al. (2008) apply a BBH model to tiie data for S5 
1803H-784 presented and discussed in this paper. The wiggling 
of the jet is attributed to a precession of the accretion disk, the 
orbital motion of the binary, and the movement of the whole 
system around the center of gravity of the AGN host (the 
"two-fluid" model is used again). 



6.1.3. Garden hose models 

Another group of models, usually identified as "garden-hose" 
models, have been extensively used for stellar jets. They assume 
a ballistic evolution of the jet, and place specific conditions 
on the distribution of the density of the jet and surrounding 
medium. Most prominent examples for this type of model are 
Hjellming & Johnston (1981) (a kinematic model for SS433 
X-ray binary), Raga et al. (1993) (an analytical model for the 
expansion of jets that get ejected with a time-dependent direc- 
tion). Biro et al. (1995) (a 2D numerical simulation based on 



the analytical model of Raga et al. (1993), considering adiabatic 
and non-adiabatic expansion of the jet), Cliffe et al. (1996) (a 
3D numerical simulation for precessing stellar jets considering 
interaction with the surrounding medium), and Stirling et al. 
(2002). Along the lines of Hjellming & Johnston (1981), Gower 
et al. (1982) produced a similar model for radio AGN jets, 
taking into account projection, relativistic, and cosmological 
redshift effects. This family of models again does not answer 
the problem of how the original precession arises, but they do 
offer a plausible explanation of the evolution of such perturbed 
jets. As all these models assume ballistic behavior, it is possible 
that they might not describe AGN jets sufficiently weU. 

Rieger (2004) presented a comparison of different configu- 
rations for curved jets, explained by balhstic and non-balhstic 
models. He proposes that a non-ballistic description of AGN jets 
is more appropriate. Geodesic (Lens-Thirring) precession of the 
accretion disk (Begelman et al. 1980) would produce a period 
much longer than observed in the light-curves of AGN (e.g. OJ 
287, S5 1803H-784). Instead those authors favor non-relativistic 
precession which leads to significantly smaller periods. Finally, 
the complete absence of magnetic fields and their effects on the 
evolution to the jet is another reason to doubt the applicability 
of the above "garden-hose" models in their present form to 
modeling of AGN jets. A model that would take into account 
MHD instabilities and the role of large-scale magnetic fields on 
the collimation and evolution of AGN jets could provide a good 
explanation of the observed structure in S5 1803H-784. We are 
currently in the process of preparing a model in this direction 
and will apply it to the data for S5 1803-1-784 presented in this 
paper. This model shall explain the position angle changes as 
well as the secondary "osciUatory" motion of the individual jet 
components. 



6.1 .4. External causes of jet curvature 

FinaUy we mention another group of models that rely on 
external causes for explanation of the curving and wiggling of 
jets. The bending of the jet could be attributed to interaction 
with the surrounding medium. This however can not fully 
explain the case of S5 1803H-784, since such periodic and fine 
structure within the jet would imply an equally finely structured 
surrounding medium, for which there is no evidence. Blandford 
& Icke (1978) attributed the structure of the jet of 3C 31 to 
dynamical interaction with a nearby galaxy (again treating the 
jet baUistically). Jagers & de Grijp (1985), based on the model 
of Blandford & Icke (1978), proposed a similar model for 
3C 10, while studying the movement of discrete blobs along 
the jet. Lupton & Gott (1984), proposed thai the curving of 
the jet of 3C 449 is due to the orbital motion of the the central 
engine around the host's center of gravity. But Ukewise, the 
periodic and finely structured curvature of the S5 1803H-784 jet 
would seem to preclude such longer-term motions as the full 
explanation. 

Yokosawa & Inoue (1985) combine some of the above 
models to the case of 3C 75, which has a double radio source 
with two pair of jets that seem to intertwine. Again treating 
the jets baUistically, the authors take into account the pressure 
of the surrounding medium, assume a time-dependent ejection 
velocity for the jets, and account for the orbital motion of the 
central engines of the jets. 
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A viable model should offer a self-consistent description and 
physical basis for the observed structure and behavior. It is im- 
portant to include the accretion disk as well as the jet into the 
model. Moreover, in the case of S5 1803H-784, the observed be- 
havior in which most jet components remain at roughly constant 
core separations should also be explained. In this light, a twist- 
ing jet model can connect the observed structure and evolution of 
the jet to the base of the jet and the accretion disk itself. We fol- 
low Steffen et al. (1995) in describing an expanding synchrotron 
blob ejected from a rotating hot spot in an accretion disk around 
a supermassive black hole. This model has been introduced to 
capture the kinematic aspects of the jet model of Camenzind & 
Krockenberger (1992) in the case of the apparent superluminal 
motion of jet components in 3C345. In our scenario the jet is 
generated by coUimated, relativistic outflow from an accretion 
disk, which is assumed to be in Keplerian rotation. Although our 
knowledge of the nozzle region between the disk and jet is not at 
all concrete, we can nonetheless make some viable assumptions 
that lead to a fairly good fit of the model to the data available. 
We are currently in the process of applying a twisted jet model 
to the data presented in this paper. Further studies on how such a 
hot spot in the accretion disk could come into being are however 
needed. 



6.2. Jet components at constant core separation 

As illustrated by the kinematic results for S5 1803H-784 pre- 
sented here, components that show no significant apparent 
superluminal (outward) motion over long time-scales can exist 
along the jet. Such components are also observed in a growing 
number of AGN (see e.g., Jorstad et al. 2001; Gomez et al.; 
Britzen et al. 2007). Different physical models and scenarios 
have been invoked to explain this phenomenon. A prominent 
example is the superluminal radio source 4C39.25, in which a 
stationary component was observed by Alberdi et al. (2000). 
They explain this in geometrical terms, proposing that the 
stationary component is at a 90° bend in the jet away from 
the observer. They disfavor the case of a recoUimation shock 
based on polarization data. Gomez et al. (2001) find three quasi- 
stationary components in the jet of 0735H-178 in a quiescent state 
of the source, and attribute them to recoUimation shocks and 
differential Doppler boosting. Britzen et al. (2009) show that 
0735-1-178 is a "hybrid" in the sense that the kinematic properties 
change simultaneously with morphological changes. At times 
when the pc-scale jet is straight, the jet components appear to be 
stationary. At times when the pc-scale jet appears to be strongly 
bent, the jet components show apparent superluminal motion. 
Britzen et al. (2007) and Kellermann et al. (2004) find stationary 
components in their 6 and 2cm AGN surveys. We also note 
the case of the quasar 3C 395 (Simon et al. 1988; Lara et al. 
1994; Lara et al. 1999), whose jet exhibits a possible stationary 
component due to possible bending of the jet. Marti & Miiller 
(2003) model AGN jets and find that pressure mismatches exist 
between the jet and the surrounding medium. These lead to 
the production of reconfinement shocks and energy density 
enhancements downstream of these shocks, which in turn give 
rise to stationary radio knots. Numerical simulations by Agudo 
et al. (2001) of an expanding relativistic jet show that the 
ejection of a superluminal component into the jet can trigger 
instabilities that can give rise to stationary components (shocks) 
in the wake of the superluminal component, especially in the 
near vicinity of the core. 



In the case of S5 1803H-784 we observe components that 
show no long-term change in their core separations, but do have 
position angles that vary in what appears to be a (quasi-)periodic 
fashion. Correlations among circular Gaussian parameters such 
as flux density, position angle, and core separation argue in favor 
of a geometric origin for the observed phenomena. We conclude 
that the kinematics in S5 1803H-784 can not be explained easily 
with recourse to published models. With curved jets, there will 
be some cases where the jet crosses the observer's line of sight, 
creating inward and stationary features. However, it is difficult 
to apply such a scenario here, in which a number of compo- 
nents show roughly constant core separations and varying posi- 
tion angles. As mentioned in the previous sub-section, we are 
currently working on a jet model to explain these characteristics 
of S5 1803H-784. 
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Clean RR map. Array; BFHKLMNOPS 
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1803 + 784 at 15.331 GHz 1996 May 17 
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Clean map. Array: BFHKLMNOPS 
1803+784 at 15,335 GHz 1997 Mar 13 
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Clean LL map. Array: BFHKLMNOPS 
1803 + ^84 at 15,335 GHz 1098 Nov 01 
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Clean map. Array: BFHKLMNOPS 
1803 + 784 at 15.335 GHz 1998 Mor 20 
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Clean map. Array: BFHKLMNOPS 
1803 + 784 at 15.363 GHz 1999 Jul 27 
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Fig. 15. Results of hybrid imaging (contours) and circular Gaussian model fitting (circles with crosses) for S5 1 803+784 at 1 5 GHz. 
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Clean LL map. Array; BFHKLMNOPS 
1803+784 at 15.335 GHz 1999 Nov 06 
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Clean LL map. Array: BFHKLMNOPS 
1803+784 at 15.335 GHz 2002 Feb 18 
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Clean map. Array: BFHKLIvlNOPS 
1803 + 784 at 15.353 GHz 2000 Jun 15 
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Clean I map. Array: BFHKLMNOPS 
1803 + 784 at 15.365 GHz 2003 Feb 06 
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Clean 1 map. Array: BFHKLMNOPS 
1803+784 at 15.365 GHz 2004 Apr 10 
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Clean 1 map. Array; BFKLNOPHM 
1803 + 784 at 15.365 GHz 2005 Apr 21 
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Clean 1 map. Array: BFHKLMNOPS 
1803+784 at 15.365 GHz 2005 Sep 05 
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(g) 



Fig. 16. Results of hybrid imaging (contours) and circular Gaussian model fitting (circles with crosses) for S5 1803+784 at 15 GHz; 
continuation of figure [TS] 
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Clean map. Array: BFHNLOSLY 
1803 + 784 at 8.413 GHz 1993 Nov 19 
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Clean map. Array: BFHKLMNOPS 
1803+784 at 8.419 GHz 1997 Dec 06 
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Clean map. Array: BFKLMNOPY 
1803 + 784 at 8.433 GHz 1997 Nov 23 




5 

Right Ascension (mas) 
Mop center: RA: 18 00 45.684, Dec: +78 28 04.018 (2000.0) 

Map peak: 2 34 Jy/beam 

Contours %: -1." 1.^ zT'T. :..S3 5.74 8.61 12.9 19.4 

Contours % 2-Z" 43 55 4 

Beom FWHM: 3.68 x 0.959 (mas) at 41,2° 



(b) 



Clean RE map. Array: BFHKLMNOPS 
1803 + 784 at 8.419 GHz 1999 May 28 
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Clean ER map. Array: BFHKLMNOPS 
1803+^84 at 8.419 GHz 2001 Feb 04 
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Fig. 17. Results of hybrid imaging (contours) and circular Gaussian model fitting (circles with crosses) for S5 1803+784 at 8.4 GHz. 
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Clean LL map. Array: BFHKLMNN0P5Y 
1803 + :"84 at 4.960 GHz 1996 Jun IT" 




Clean LL map. Array: BEFHKLMNOPSY 
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Clean LL map. Array: BEFHKKLMNOPS 
1803 + :"64 at 4.800 GHz 1999 May OS 




Clean map. Array: BEFHJKLMNNOPSVVY 
1803 + 784 at 4,965 GHz 1999 Jun 10 
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Clean LL map. Array: BFHKLNOPRSWYM 
J1800+78 at 4.966 GHz 2001 Apr 17 
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Fig. 18. Results of hybrid imaging (contours) and circular Gaussian model fitting (circles with crosses) for S5 1803+784 at 5 GHz. 
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Clean LL map. Array; DFHNMNOPBKMSL 
1803+784 at 2.287 GHz 1995 Jun 03 
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Fig. 19. Results of hybrid imaging (contours) and circular Gaussian model fitting (circles with crosses) for S5 1803+784 at 2.3 GHz. 




Fig. 20. Results of hybrid imaging (contours) and circular Gaussian model fitting (circles with crosses) for S5 1 803+784 at 1 .6 GHz. 
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Table 1. Circular Gaussian model-fit results for S5 1803+784. The columns list the epoch of observation, the jet-component identifi- 
cation, the flux density, the radial distance of the component center from the core, the position angle of the center of the component, 
the FWHM major axis, and the chi-squared of the fit. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Epoch 


Id. 


S[Jy] 


r [mas] 


e [deg] 


M.A. [mas] 




1.6 GHz 


1998.92 


r 


0.946+0.142 


0.00 


0.0 


0.26+0.02 


38.232803 


1998.92 


Ca 


0.558+0.084 


1.41+0.05 


262.7+0.8 


0.69+0.03 




1998.92 


C4 


0.091+0.014 


4.35+0.44 


267.8+1.0 


2.82+0.30 




1998.92 


C8 


0.059±0.009 


7.99+0.80 


265.4+2.0 


1.30+0.35 




1998.92 


C30 


0.259+0.052 


29.67+3.00 


256.7+4.0 


15.72+4.00 




2001.87 


r 


1.233+0.169 


0.00+0.01 


0.0+0.5 


0.60+0.08 


0.0923 


2001.87 


Ca 


0.266+0.047 


2.33+0.71 


-94.8+1.5 


1.12+0.29 




2001.87 


C8 


0.134+0.019 


6.15+0.98 


-96.7+1.2 


3.17+0.57 




2001.87 


C12 


0.058+0.010 


9.45+0.42 


-97.0+0.9 


3.70+0.11 




2001.87 




0.010+0.001 


16.85+1.10 


-91.4+0.9 


5.10+1.60 




2001.87 


C30 


0.241+0.034 


28.94+0.49 


-104.7+1.0 


14.18+0.79 




2001.87 




0.045+0.007 


43.07+0.61 


-100.0+0.8 


15.98+3.20 




2.3 GHz 


1996.42 


r 


1.642+0.246 


0.00 


0.0 


0.00+0.01 


86.988178 


1996.42 


Ca 


0.481+0.072 


1.43+0.07 


-92.2+0.5 


0.01+0.01 




1996.42 


C4 


0.041+0.006 


3.91+0.39 


-85.3+1.0 


0.44+0.02 




1996.42 


C8 


0.124+0.025 


6.71+0.80 


-89.4+1.8 


3.19+0.32 




1996.42 


C12 


0.028+0.010 


11.92+2.38 


-102.7+1.9 


1.47+0.15 




1996.42 


C3() 


0.204+0.069 


26.67+3.00 


-102.2+3.9 


13.47+2.70 




5.0 GHz 


1996.46 


r 


1.673+0.251 


0.00 


0.0 


0.36+0.02 


52.027293 


1996.46 


Ca 


0.422+0.063 


1.35+0.07 


-96.8+1.1 


0.86+0.04 




1996.46 




0.090+0.014 


5.58+0.60 


-89.9+1.54 


4.44+0.44 




1997.15 


r 


1.302+0.195 


0.00 


0.0 


0.01+0.01 


27.458208 


1997.15 


CI 


0.274+0.041 


0.87+0.04 


-94.2+0.9 


0.33+0.07 




1996.15 


Ca 


0.218+0.033 


1.73+0.09 


-97.3+0.7 


0.53+0.08 




1997.15 


C4 


0.017+0.002 


3.95+0.40 


-85.0+1.4 


0.47+0.10 




1997.15 


C8 


0.081+0.016 


6.36+0.66 


-99.7+1.9 


4.42+0.44 




1997.15 


C30 


0.067+0.013 


27.87+2.81 


-107.6+3.9 


6.48+0.84 




1999.35 


r 


1.774+0.266 


0.00 


0.00 


0.30+0.02 


15.738181 


1999.35 


CI 


0.033+0.005 


0.77+0.04 


-69.4+0.5 


0.00+ 





1999.35 Ca 
1999.35 C4 



0.507+0.189 
0.031+0.025 



1.31+0.07 
3.46+0.35 



-93.0+1.5 
-98.1 + 1.1 



0.81+0.12 
0.81+0.41 



1999.44 r 
1999.44 CI 
1999.44 Ca 
1999.44 C2 
1999.44 C4 
1999.44 C12 
1999.44 C30 



1.423+0.213 
0.450+0.067 
0.274+0.041 
0.101+0.015 
0.071+0.011 
0.033+0.006 
0.099+0.020 



0.00 
0.53+0.03 

1.34+0.07 
1.86+0.09 
4.81+0.48 

10.00+1.00 
28.27+2.83 



0.0 

-75.7+1.0 
-89.2+0.7 
■100.6+0.8 
-94.8+1.1 

-93.9+1.6 
■104.2+2.7 



0.21+0.01 
0.22+0.01 

0.49+0.03 
0.78+0.08 
2.94+0.30 

1.99+0.20 
14.06+2.58 



2001.29 r 
2001.29 CI 
2001.29 Ca 
2001.29 C4 
2001.29 C8 
2001.29 C12 
2001.29 C30 



1.602+0.246 
0.511+0.071 
0.228+0.069 
0.045+0.023 
0.052+0.022 
0.039+0.018 
0.112+0.018 



0.00+0.01 
0.85+0.08 
1.62+0.38 
3.89+0.18 
6.07+0.24 
9.23+0.35 
29.12+0.59 



0.0+0.5 
-81.1+2.1 
-91.7+3.5 
-84.7+4.9 
-98.3+0.8 
-95.1 + 1.7 
■103.3+1.2 



0.21+0.05 
0.43+0.11 
0.72+0.22 
1.38+0.61 
2.05+0.85 
3.22+0.68 
11.83+0.19 



0.0477 



8.4 GHz 


1993.88 


r 


1.390+0.208 


0.00 


0.00 


0.01+ 2.8677791 


1993.88 


CO 


0.014+0.002 


0.25+0.01 


-69.3+0.7 


0.00+0.01 


1993.88 


Cx 


0.262+0.039 


0.33+0.02 


-99.5+0.4 


0.00+0.01 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Epoch 


Id. 


S[Jy] 


r [mas] 


6 [deg] M.A. [mas] 


T 


1993.88 


Ca 


0.263+0.039 


1.45+0.07 


-91.8+0.8 


0.36+0.02 




1993.88 


C4 


0.043±0.016 


3.07+0.11 


-87.5+1.2 


1.03+0.17 




1997.90 


r 


2 467±0 370 


0.00 


0.00 


0.45+0.02 


6 8563583 


1997.90 


Ca 


719±0 108 


1 39±0 07 


-97.4+1.3 


0.68+0.03 




1997.90 


C4 


127+0 019 


4 00+0 45 


-89.5+1.3 


1.31+0.13 




1997.93 


I* 


1 ^-s^+o 203 


0.00 


0.00 


0.12+0.01 


1 7453468 


1997.93 


CO 


270+0 041 


28+0 01 


-83.5+1.9 


0.36+0.02 




1997.93 


CI 


0.176±0.026 


1.01+0.05 


-86.6+0.6 


0.44+0.02 




1997.93 


Ca 


0.212+0.032 


1.54+0.08 


-98.0+0.5 


0.45+0.02 




1997.93 


C2 


0.035+0.005 


2.23+0.11 


-97.0+1.0 


0.76+0.04 




1997.93 


C4 


0.004+0.001 


3.16+0.32 


-109.5+1.0 


0.28+0.03 




1997.93 


C4 


0.028+0.004 


3.89+0.39 


-89.5+1.0 


1.38+0.14 




1997.93 


C8 


0.039+0.008 


7.71+0.90 


-96.6+3.0 


3.00+0.80 





1999.41 
1999.41 
1999.41 
1999.41 
1999.41 
1999.41 
1999.41 
1999.41 
1999.41 



r 1, 
CO 
CI 
Ca 
C2 
C4 0, 
C8 
C12 
C30 

r 1, 
CO 
CI 
Ca 
C2 
C4 
C8 
C12 
C30 



257+0.188 
,496+0.074 
,171+0.026 
.140+0.021 
086+0.013 
046+0.007 
025+0.004 
.020+0.003 
036+0.006 



0.00 
0.33+0.02 
0.92+0.05 
1.50+0.08 
1.84+0.09 
3.43+0.34 
6.47+0.65 
9.40+0.90 
27.91+2.40 



0.0 

-79.7+1.1 
-78.9+0.5 
-92.5+0.5 
-96.9+0.5 
-93.9+1.0 
-95.7+0.9 
-95.6+1.0 
■103.1+3.2 



0.12+ 
0.28+0.01 
0.43+0.02 
0.32+0.02 
0.71+0.04 
1.74+0.09 
2.43+0.24 
2.96+0.30 
7.51+0.75 



1.5628510 



2001.09 
2001.09 
2001.09 
2001.09 
2001.09 
2001.09 
2001.09 
2001.09 
2001.09 



317+0.146 
437+0.068 
.259+0.094 
.179+0.032 
059+0.008 
055+0.010 
027+0.003 
024+0.008 
.018+0.003 



0.00+0.01 
0.37+0.07 
0.93+0.12 
1.47+0.04 
1.95+0.04 
4.11+0.33 
6.12+0.15 
8.95+0.73 
27.61+0.12 



0.0+0.5 
-81.0+1.1 
-79.4+1.7 
-91.6+0.8 
-93.7+0.6 
-91.0+0.9 
-96.5+0.3 
-94.1+0.5 
-103.8+0.3 



0.13+0.09 
0.30+0.03 
0.42+0.19 
0.45+0.11 
0.72+0.04 
2.44+0.64 
2.29+0.13 
2.90+0.37 
7.18+0.19 



1.1756 



15 GHz 



1994.67 


r 


1.115+0.167 


0.00 


0.00 


0.13+0.01 0.43219569 


1994.67 


CO 


0.176+0.030 


0.39+0.03 


-92.8+3.9 


0.40+0.05 


1994.67 


Ca 


0.097+0.015 


1.46+0.04 


-93.9+1.5 


0.61+0.05 


1994.67 


C4 


0.030+0.010 


3.43+0.10 


-90.6+1.7 


1.57+0.12 



1996.38 


r 


1.559+0.234 


0.00 


0.00 


0.12+0.01 0.43219569 


1996.38 


CO 


0.572+0.086 


0.25+0.01 


-93.7+2.1 


0.23+0.01 


1996.38 


CI 


0.094+0.014 


0.67+0.03 


-98.0+1.1 


0.35+0.02 


1996.38 


Ca 


0.191+0.029 


1.44+0.07 


-92.9+0.6 


0.45+0.02 


1996.38 


C2 


0.124+0.019 


1.74+0.09 


-96.1±0.5 


0.45 ±0.02 


1996.38 


C4 


0.022+0.003 


3.38+0.34 


-86.6+1.0 


1.43+0.14 



1997.20 


r 


1.362+0.204 


0.00 


0.00 


0.11+0.01 1.4601196 


1997.20 


CO 


0.295+0.044 


0.30+0.02 


-89.0+1.6 


0.31+0.02 


1997.20 


CI 


0.135+0.020 


0.84+0.04 


-90.3+0.9 


0.55+0.03 


1997.20 


Ca 


0.111+0.017 


1.47+0.07 


-96.5+0.5 


0.39+0.02 


1997.20 


C2 


0.101+0.015 


1.78+0.09 


-96.3+0.5 


0.59+0.03 


1997.20 


C4 


0.030+0.005 


3.13+0.31 


-90.9+1.0 


1.64+0.16 


1997.20 


C8 


0.034+0.005 


6.81+0.68 


-97.5+1.0 


3.40+0.34 


1998.22 


r 


1.089+0.163 


0.00 


0.00 


0.25+0.01 0.44147467 


1998.22 


CO 


0.916+0.137 


0.23+0.01 


-77.2+5.3 


0.30+0.02 


1998.22 


CI 


0.119+0.018 


0.93+0.05 


-83.7+1.2 


0.47+0.02 


1998.22 


Ca 


0.247+0.037 


1.49+0.07 


-95.0+0.7 


0.46+0.02 


1998.22 


C2 


0.049+0.007 


2.04+0.20 


-96.8+1.0 


0.62+0.06 


1998.22 


C4 


0.038+0.008 


3.40+0.34 


-89.6+1.4 


2.07+0.21 


continued on next page 



S. Britzen et al.: The kinematics in the pc-scale jets of AGN 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Epoch 


Id. 


S[Jy] 


r [mas] 


[deg] 


M.A. [mas] 


T 


1998.84 


I* 


1.113±0.167 


0.00 


0.00 


08+0 01 


1 4787256 


1998.84 


CO 


650+0 098 


18+0 02 


-86.2+9.4 


35+0 02 




1998.84 


CI 


102±0 015 


92+0 05 


-75.6+2.8 


60+0 03 




1998.84 


Ca 


0.141±0.021 


1.42+0.07 


-92.4+0.6 


0.22+0.01 




1998.84 


C2 


0.147±0.022 


1 67+0 08 


-99.0+1.3 


53+0 03 




1999.57 


I* 


2 01 5+0 102 


0.00 


0.0 


1 2+0 01 


1 1925775 

-L.J- J7^^ 1 1 ,J 


1999.57 


CO 


430+0 064 


33+0 02 


-77.0+2.1 


26+0 01 

\J , i-i\J -!—\J ,\J X 




1999.57 


CI 


1 78+0 027 


84+0 04 


-77.1+0.9 


31+0 02 




1999.57 


Ca 


048+0 007 


1 ^2+0 07 

J. . ^ ^ -1—\J . \J 1 


-83 5+0 6 


0.01 + 




1999.57 


C2 


0.131+0.020 


1.57+0.08 


-92.5+0.6 


0.42+0.07 




1999.57 




040+0 006 


1 98+0 10 


-93.4+1.0 


68+0 11 




1999.57 


C4 


032±0 005 


3 28+0 33 


-92.8+1.1 


1.57+0.14 




1999.57 




007+0 002 


5.14+0.51 


-88 7+2 


1.01+0.18 




1999.57 


C8 


020+0 004 


7.38+0.74 


-96.2+2.0 


2.47+0.34 




1999.57 


C12 


006+0 001 


10.26+1.16 


-91.7+2.0 


3.04+0.71 




1999.85 


I* 


1 603+0 241 

X .\J\J ^ _!_ \J . J. 


0.00 


0.0 


090+0 01 


1 1036355 


1999.85 


CO 


449+0 067 


28+0 05 


-79.6+2.9 


24+0 01 




1999.85 


CI 


0.168+0.024 


0.94+0.06 


-76.2+1.8 


0.39+0.02 




1999.85 


Ca 


093+0 015 


1.38+0.17 


-86 8+0 6 


0.17+0.01 




1999.85 


C2 


098+0 01 5 

\J t\J ^ yJ -!—\J ,\J J. ^ 


1 64+0 10 


-94.5+0.5 


42+0 02 




1999.85 


C4 


039±0 007 


2.44+0.13 


-97.5+2.8 


1.27+0.13 




1999.85 




0.021+0.004 


4.55+0.12 


-88.1+1.0 


1.61+0.16 




1999.85 


C8 


013+0 003 


7 86+0 1 8 


-97.6+2.0 


2 05+0 21 




2000.46 


r 


1.606+0.241 


0.00 


0.00 


0.13+0.01 


0.46644555 


2000.46 


CO 


0.364+0.055 


0.35+0.02 


-79.2+1.8 


0.27+0.01 




2000.46 


CI 


0.144+0.022 


0.96+0.05 


-79.3+0.7 


0.36+0.02 




2000.46 


Ca 


0.134+0.020 


1.44+0.07 


-88.8+0.5 


0.29+0.02 




2000.46 


C2 


0.105+0.016 


1.70+0.03 


-92.1+0.7 


0.62+0.03 




2000.46 


C4 


0.055+0.008 


3.61+0.72 


-91.9+1.1 


1.95+0.20 





1.3686 



2002.13 r 1.702+0.251 0.00 + 0.01 0.0+0.3 0.11+0.02 

2002.13 CO 0.151+0.023 0.37+0.01 -94.6+1.2 0.29+0.01 

2002.13 B3 0.008+0.001 0.76+0.02 -89.8+ 1.2 0.14+0.03 

2002.13 CI 0.072+0.012 1.02+0.02 -75.2+0.5 0.35+0.02 

2002.13 Ca 0.318+0.043 1.43+0.02 -88.1+0.5 0.48+0.05 

2002.13 C2 0.029+0.004 2.02+0.02 -94.0+0.5 0.68+0.01 

2002.13 C4 0.031+0.005 3.61+0.05 -89.0+0.8 1.94+0.02 

2002.13 C8 0.044+0.007 6.37+0.09 -95.1+0.8 3.40+0.02 

2002.13 C12 0.009+0.001 11.36+0.37 -94.6+ 1.9 3.37+0.67 



2003. 


10 


r 


1.245+0.179 


0.00+ 0.02 


0.0+ 0.2 


0.10+0.03 1.5375 


2003. 


10 


CO 


0.310+ 0.048 


0.31+0.03 


-76.9+ 1.7 


0.27+0.01 


2003. 


10 


CI 


0.039+ 0.008 


0.66+ 0.03 


-84.2+ 1.3 


0.09+0.04 


2003. 


10 


Ca 


0.287+ 0.058 


1.42+0.01 


-88.4+ 0.3 


0.35+0.09 


2003. 


,10 


C2 


0.113+0.045 


1.78+0.24 


-89.1+0.8 


1.00+0.26 


2003. 


10 


B3 


0.035+0.016 


2.64+ 0.38 


-88.9+ 1.2 


1.17+0.19 


2003. 


10 


C4 


0.039+ 0.006 


4.67+ 0.32 


-90.0+ 0.8 


2.01+0.35 


2003. 


10 


C8 


0.031+0.004 


7.44+ 0.49 


-97.1+0.7 


3.00+0.16 


2003. 


,10 


C12 


0.015+ 0.004 


12.37+0.78 


-93.2+0.6 


4.31+0.64 



1.4844 



2004.27 r 1.264+0.185 0.00+0.01 0.0+0.1 0.09+0.03 

2004.27 CO 0.117+0.017 0.33+0.01 -89.1+ 1.3 0.25+0.01 

2004.27 CI 0.124+0.018 1.00+0.05 -78.0+0.6 0.38+0.05 

2004.27 Ca 0.185+0.028 1.47+0.02 -90.6+0.4 0.31+0.06 

2004.27 C2 0.046+0.009 1.76+0.06 -95.6+0.5 0.49+0.04 

2004.27 C4 0.044+0.007 3.64+0.09 -86.9+0.9 1.83+0.05 

2004.27 C8 0.046+0.007 6.42+0.24 -96.3+0.8 2.71+0.08 

2004.27 C12 0.014+0.003 9.81+0.30 -96.3+0.6 3.09+0.55 
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(1) (2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Epoch Id. 


S[Jy] 


r [mas] 


[deg] 


M.A. [mas] 





2005.30 r 0.819+0.119 0.00+0.01 0.0+0.2 0.12+0.04 0.9761 

2005.30 CO 0.160+0.028 0.22+0.02 -61.6±0.3 0.05+0.03 

2005.30 CI 0.091± 0.014 0.61±0.02 -75.8± 1.3 0.35±0.02 

2005.30 Ca 0.297± 0.045 1.47±0.01 -87.4±0.6 0.42±0.03 

2005.30 C2 0.040+0.011 1.98±0.06 -84.1+0.5 0.51±0.07 

2005.30 C4 0.018+0.005 3.98+0.09 -84.9+0.8 1.12+0.34 

2005.30 C8 0.064+0.009 6.08+0.30 -94.7± 1.0 3.04+0.23 

2005.30 €12 0.010+0.005 9.80+0.97 -91.9+ 2.5 2.90+0.70 



2005 


.68 


r 


1.381+0.203 


0.00+ 0.01 


0.0+ 0.3 


0.09+0.02 0.8685 


2005 


.68 


CO 


0.146+0.022 


0.25+ 0.03 


-74.2+ 1.5 


0.21+0.02 


2005 


.68 


CI 


0.048+ 0.008 


0.75+ 0.02 


-78.0+ 0.8 


0.32+0.03 


2005 


.68 


Ca 


0.244± 0.064 


1.44± 0.02 


-89.0± 0.3 


0.31 ±0.09 


2005 


.68 


C2 


0.120+ 0.019 


1.71+0.04 


-88.8+ 0.5 


0.55+0.02 


2005 


.68 


C4 


0.021+0.007 


3.44+ 0.29 


-87.0+ 1.3 


1.50+0.20 


2005 


.68 


B3 


0.013+0.001 


4.72+ 0.06 


-83.7+0.4 


1.16+0.03 


2005 


.68 


C8 


0.058+ 0.006 


6.37+0.11 


-97.4+ 0.6 


2.49+0.18 


2005 


.68 


C12 


0.022± 0.001 


8.98± 0.20 


-93.7± 0.7 


3.68±0.22 
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Table 2. Gaussian model-fit results for S5 1803H-784 obtained from the literature. The columns list the epoch of observation, the 
jet-component identification, the flux density, the radial distance of the component center from the core, the position angle of the 
center of the component, the FWHM major axis, and for elliptical Gaussian fits, the axial ratio and the position angle of the major 
axis. The final column provides a reference to the original pubUcation. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Epoch 


Id. 


S[Jy] 


r [mas] 


e [deg] 


Ma.A. [mas] 


Axial ratio 


<I> [deg] 


Reference 


2.32 GHz 


1994.52 


1 


1.63 


0.00 


0.00 


1.37 


0.00 


87 


Fey et al. 1996 


1994.52 


2 


0.34 


4.0 


-94 


7.24 


0.28 


81 




1994.52 


3 


0.28 


27.6 


-105 


12.11 


1.00 






5 GHz 


1997.79 


CI 


0.413 


1.38 


-98 


0.65 


0.82 


50 


Lister et al. 2001 


8 GHz 


1988.79 


Core 


2.45 


0.00 


0.0 


0.1 






Tateyama et al. 2002 


1988.79 


CI 


0.23 


1.50 


-95 


0.1 








1988.79 


C2 


0.21 


0.90 


-95 


0.1 








1989.96 


Core 


2.77 


0.00 


0.0 


0.1 






Tateyama et al. 2002 


1989.96 


CI 


0.22 


1.51 


-94 


0.1 








1989.96 


C2 


0.20 


1.09 


-90 


0.3 








1989.96 


E 


0.65 


0.40 


-81 


0.2 








1991.96 


Core 


1.38 


0.00 


0.0 


0.1 






Tateyama et al. 2002 


1991.96 


CI 


0.14 


1.71 


-90.9 


0.1 








1991.96 


C2 


0.27 


1.28 


-87.4 


0.1 








1991.96 


E 


0.06 


0.58 


-85.1 


0.1 








1992.62 


Core 


1.68 


0.00 


0.0 


0.2 






Tateyama et al. 2002 


1992.62 


CI 


0.11 


1.80 


-90 


0.1 








1992.62 


C2 


0.40 


1.27 


-95 


0.1 








1992.62 


E 


0.32 


0.48 


-93 


0.01 








1993.87 


Core 


1.25 


0.00 


0.0 


0.2 






Tateyama et al. 2002 


1993.87 


C2 


0.33 


1.39 


-97.2 


0.1 








1994.71 


Core 


1.49 


0.00 


0.0 


0.2 






Tateyama et al. 2002 


1994.71 


C2 


0.17 


1.43 


-96.2 


0.1 








1995.71 


Core 


2.10 


0.00 


0.0 


0.25 






Tateyama et al. 2002 


1995.71 


C2 


0.28 


1.52 


-93.0 


0.1 








1995.71 


C3 


0.18 


1.12 


-93.1 


0.1 








1996.29 


Core 


1.57 


0.00 


0.0 


0.2 






Tateyama et al. 2002 


1996.29 


C2 


0.10 


1.59 


-91 


0.1 








1996.29 


C3 


0.04 


1.10 


-97 


0.1 








1996.29 


E 


0.18 


0.44 


-101 


0.1 








1996.71 


Core 


1.56 


0.00 


0.0 


0.2 






Tateyama et al. 2002 


1996.71 


C2 


0.15 


1.73 


-92.4 


0.1 








1996.71 


C3 


0.10 


1.25 


-96.1 


0.1 








1996.71 


E 


0.17 


0.51 


-89.2 


0.1 








1998.79 


Core 


1.60 


0.00 


0.0 


0.2 






Tateyama et al. 2002 


1998.79 


C2 


0.13 


1.82 


-97.0 


0.1 








1998.79 


C3 


0.26 


1.30 


-93.7 


0.1 








1998.79 


E 


0.41 


0.40 


-93.3 


0.1 








1999.62 


Core 


1.89 


0.00 


0.0 


0.2 






Tateyama et al. 2002 


1999.62 


C2 


0.09 


1.86 


-90.0 


0.1 








1999.62 


C3 


0.19 


1.33 


-88.1 


0.1 








1999.62 


E 


0.28 


0.62 


-77.8 


0.1 








8.4 GHz 


1997.93 


XA 


1.493±0.001 


0.00 




0.2±0.1 


0.8+0.1 


+119+1 


Rosetal. 2001 


1997.93 


XB 


0.231±0.001 


0.6±0.2 


-86±19 


0.7±0.1 


0.6+0.1 


+102+1 




1997.93 


XC 


0.274±0.001 


1.4±0.2 


-97±8 


0.6±0.1 


0.7±0.1 


+41±1 




1997.93 


XD 


0.036+0.001 


2.0+0.2 


-98+5 


0.5+0.1 


0.2+0.1 


-12+4 




1997.93 


XE 


0.024+0.001 


2.7+0.2 


-97+4 


1.4+0.1 


0.4+0.1 


+ 1+3 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Epoch 


Id. 


S[Jy] 


r [mas] 


e [deg] 


Ma.A. [mas] 


Axial ratio 


O [deg] 


Reference 


1997.93 


XF 


018+0 001 


3.9±0.2 


-89±3 


1.1 ±0.1 


0.9+0.1 


-30+16 




1997.93 


XG 


019+0 001 


6.3+0.2 


-96+3 


2.8+0.2 


0.5+0.1 


+31+4 




1997.93 


XH 


025+0 001 


8.9+0.2 


-95 ±1 


4.1+0.2 


0.5+0.1 


+99+4 




1999.41 


XA 


1 517+0 001 

t- ,,J t- 1 J^\J ,\J\J J- 


0.00 




0.2±0.1 


0.6±0.1 


+107±1 


Ros et al 2001 


1999.41 


XB 


n ^74+0 001 

Vy.J / ^ J—\J .\J\J 1 


0.5+0.2 


-81+21 


0.5+0.1 


0.7+0.1 


+117+1 




1999.41 


XC 


0.192+0.001 


1.4±0.2 


-93±8 


0.5±0.1 


0.3+0.1 


+22+1 




1999.41 


XD 


0.058±0.001 


1.9+0.2 


-94+6 


0.6±0.1 


0.2+0.1 


+3+2 




1999.41 


XE 


0.030±0.001 


2.7±0.2 


-93 ±4 


1.9±0.1 


0.4±0.1 


+126±2 




1999.41 


XF 


0.030+0.001 


4.0+0.2 


-97±3 


2.4±0.1 


0.4+0.1 


+111+1 




1999.41 


XG 


0.020±0.001 


7.0±0.2 


-97±2 


2.9±0.2 


0.6±0.1 


+115+5 




1999.41 


XH 


0.017±0.001 


9.4±0.2 


-97±1 


3.6±0.3 


0.6±0.1 


+84±5 




8.55 GHz 


1994.52 


1 


1.91 


0.00 




0.19 


0.63 


-81 


Fey et al. 1996 


1994.52 


2 


0.38 


1.0 


-92 


1.53 


0.28 


83 




1994.52 


3 


0.05 


3.5 


-86 


1.26 


1.00 






1994.52 


4 


0.05 


6.6 


-95 


2.57 


1.00 






86 GHz 


1993.29 


C 


0.448+0.063 


0.00 




0.08+0.03 






Lobanov et al. 2000 


1993.29 


Jl 


0.376+0.056 


0.27+0.06 


-0.9±8.0 


0.08+0.03 
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Table 7. Discrete cross-correlations (based on components from the 8 GHz model fits) and Pearson's correlation coefficients (from 
15 GHz model fits), calculated between the time series of various pairs of jet parameters and/or the total flux densities. Column 1 
gives the frequency of the analyzed data set; colum 2 lists the individual jet component from which the parameters to correlate are 
taken; columns 3 and 4 list the parameters correlated. For 8 GHz data, the DCF peak and t are the peak of discrete cross-correlation 
function and the associated time lag in years. For 15 GHz data, Corr.coeff. is the value of the Pearson's correlation coefficient, and 
p is the probability of getting such a high correlation by chance. 



V 


Comp. 


Par. 1 


Par. 2 


DCF 


T 


Corr. 


P 


[GHz] 


ID 






peak 




coeff. 




8 


Ca 


core Sep. 


total flux 8 GHz 


-0.68 ± 0.06 


-1.12 ±0.05 










p.a. 


total flux 8 GHz 


-0.77 ± 0.04 


0.72 ± 0.06 








C2i 


core Sep. 


total flux 8 GHz 


0.56 ± 0.06 


-0.07 ± 0.08 








C23 


p.a. 


total flux 8 GHz 


0.60 ± 0.07 


-0.15 ±0.08 








CI 


core Sep. 


total flux 14.5 GHz 


0.41 ± 0.07 


-0.07 ± 0.04 










p.a. 


total flux 14.5 GHz 


0.61 ± 0.05 


0.06 ± 0.06 






8 


CI 


p.a. 


core Sep. 


0.54 ± 0.09 


0.20 ± 0.05 










p.a. 


flux 


0.53 ± 0.08 


-0.08 ± 0.04 










core Sep. 


flux 


0.41 ± 0.06 


-0.58 ± 0.06 








Ca 


p.a. 


core Sep. 


-0.59 ± 0.06 


-0.59 ± 0.05 










p.a. 


flux 


0.42 ± 0.06 


-0.56 ± 0.06 










core Sep. 


flux 


-0.64 + 0.06 


-0.03 ± 0.04 








C2i 


p.a. 


core Sep. 


-0.53 ± 0.09 


0.31 ± 0.08 








C22 


p.a. 


core Sep. 


-0.61 ± 0.08 


-0.42 ± 0.04 








C4 


p.a. 


core Sep. 


0.97 ±0.10 


0.09 ± 0.03 






15 


CO 


p.a. 


core Sep. 






-0.44 ± 0.09 


0.131 






core Sep. 


flux 






-0.52 ± 0.08 


0.067 




CI 


p.a. 


core Sep. 






0.40 ±0.10 


0.198 






core Sep. 


flux 






0.46 ± 0.09 


0.164 




Ca 


p.a. 


core Sep. 






-0.73 ± 0.05 


0.008 






core Sep. 


flux 






0.57 ± 0.07 


0.054 




C2 


core Sep. 


flux 






-0.83 ± 0.03 


0.002 




C4 


p.a. 


core Sep. 






0.58 ± 0.07 


0.064 






p.a. 


flux 






-0.51 ±0.08 


0.107 
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Table 8. Discrete cross-correlations (based on components from the 8 GHz model fits) and Pearson's correlation coefficients (from 
15 GHz model fits), calculated for various pairs jet parameters and the core flux density. Column 1 gives the frequency of the 
analyzed data set; columns 2 and 3 lists the two jet components participating in the correlation; colunrn 4 shows the circular 
Gaussian parameter correlated. When "core" is the first component, its parameter is always flux-density. For 8 GHz data, DCF peak 
and T are the peak of the discrete correlation function and the associated time lag in years. For 15 GHz data, Corr.coeff. is the value 
of the Pearson's correlation coefficient, and p the probability of getting such a high correlation by chance. 



V 


Comp. 


Comp. 


Par. 


DCF 


T 


Corr. 


P 


[GHz] 


ID 1 


ID 2 




peak 




coeff. 




8 


CO 


Ca 


core Sep. 


-0.48 ± 0.12 


0.54 ± 0.07 








CO 


C2i 


p.a. 


-0.49 + 0.14 


-0.32 ± 0.06 








CO 


C23 


p.a. 


0.69 ± 0.18 


0.14 + 0.08 








CI 


Ca 


core Sep. 


0.47 ± 0.06 


0.08 ± 0.06 












p.a. 


-0.48 ± 0.07 


0.48 ± 0.07 








CI 


C2i 


core Sep. 


-0.59 + 0.09 


-0.80 + 0.09 












p.a. 


-0.45 + 0.09 


-0.29 + 0.09 








CI 


C22 


core Sep. 


0.53 ± 0.14 


-0.02 + 0.06 












p.a. 


0.96 ± 0.03 


-0.02 + 0.02 








CI 


C23 


core Sep. 


0.83 + 0.06 


0.81 ± 0.02 












p.a. 


0.72 ± 0.09 


-0.06 ± 0.06 








Ca 


C2i 


core Sep. 


-0.48 + 0.09 


-0.09 ± 0.06 








Ca 


C23 


core Sep. 


0.62 ± 0.09 


0.55 ± 0.06 












p.a. 


-0.61 + 0.07 


0.49 ± 0.08 








C2i 


C22 


p.a. 


-0.53 + 0.08 


0.42 + 0.05 








C2i 


C23 


core Sep. 


0.62 ± 0.09 


-0.42 ± 0.06 












p.a. 


-0.78 ± 0.09 


0.17+0.02 








C22 


C23 


p.a. 


0.79 ± 0.05 


-0.21 + 0.05 






8 


Core 


CO 


core Sep. 


0.62 ± 0.09 


-0.23 ± 0.07 












p.a. 


0.77 ±0.10 


-0.09 ± 0.05 








Core 


CI 


core Sep. 


-0.47 ± 0.07 


0.35 ± 0.08 












p.a. 


-0.45 ± 0.07 


-0.31 ± 0.04 








Core 


Ca 


flux 


0.59 ± 0.06 


0.49 + 0.10 












core Sep. 


-0.91 + 0.04 


0.23 ± 0.03 












p.a. 


-0.83 + 0.04 


-0.92 ± 0.04 








Core 


C2i 


core Sep. 


0.64 ± 0.07 


0.02 ± 0.05 








Core 


C22 


p.a. 


-0.47 + 0.08 


-0.29 ± 0.05 








Core 


C2, 


core Sep. 


0.77 ± 0.05 


-0.36 ± 0.04 












p.a. 


0.54 ± 0.08 


0.37 ± 0.10 






15 


Core 


CO 


core Sep. 






0.65 ± 0.05 


0.021 




Core 


Ca 


core Sep. 






-0.72 ± 0.05 


0.008 








p.a. 






0.49 ± 0.09 


0.146 




Core 


C2 


core Sep. 






-0.49 ± 0.08 


0.124 




Core 


C4 


p.a. 






-0.56 ± 0.08 


0.074 


15 


CO 


C2 


p.a. 






0.77 ± 0.04 


0.006 




CI 


Ca 


flux 






-0.84 ± 0.03 


0.001 








p.a. 






-0.67 + 0.06 


0.025 




Ca 


C2 


flux 






-0.60 ± 0.07 


0.051 








core Sep. 






0.69 + 0.06 


0.019 








p.a. 






0.54 ± 0.08 


0.084 



